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ABSTRACT 


During the period of this research project, a comprehensive study of pyramidal horn 
antennas was conducted. Full-wave analytical and numerical techniques have been 
developed to analyze horn antennas with or without impedance surfaces. Based on 
these full- wave analytic techniques, research was conducted on the use of impedance 
surfaces on the walls of the horn antennas to control the antenna radiation patterns 
without a substantial loss of antenna gain. It was found that the use of impedance 
surfaces could modify the antenna radiation patterns. In addition to the analytical 
and numerical models, experimental models were also constructed and they were 
used to validate the predictions. Excellent agreement between theoretical predictions 
and the measured data was obtained for pyramidal horns with perfectly conducting 
surfaces. Very good comparisons between numerical and experimental models were 
also obtained for horns with impedance surfaces. 
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CHAPTER 1 


MOMENT METHOD ANALYSIS OF HORN ANTENNAS 
1.1 Introduction 

The horn antenna is the simplest and probably the most widely used microwave 
radiator. It is used as the feed for large reflectors and lens antennas in communication 
systems throughout the world. It is also a high gain element in phased arrays. Horn 
antennas are highly accurate radiating devices, and are often used as standard-gain 
devices for the calibration of other antennas. Applications of horn antennas have been 
explored for nearly a century. In addition, extensive investigations of horn antennas 
have been of increasing interest during the past three decades. Some of the early 
research papers on the horn antennas are well documented in Love’s collection[l]. 

Besides being a high-gain and high-efficiency microwave antenna, the pyramidal 
horn exhibits some additional advantages. Its rectangular geometry is easy to con- 
struct, and hence the antenna is a low-cost device. The aperture size of the horn can 
be adjusted to achieve some specific beam characteristics with negligible changes 
in other properties. Furthermore, the pyramidal horn can easily be excited with 
conventional microwave circuit devices. 

Analysis and design of the pyramided horns was traditionally conducted using 
approximate aperture field distributions. Using this approach the contributions from 
induced currents on other parts of the horn surface are often assumed negligible. 
A quadratic phase term is usuadly assumed to account for the flaring of the horn 
transition [2, 3, 4]. This approximate method predicts fairly well the main-beam 
of the far-field radiation pattern and the gain of the antenna. Since reflections, 
mode couplings, and diffracted fields from the exterior surfaces are not included, the 
approximate aperture field method does not predict very well the sidelobes in the 
back region of the antenna. 

In the 1960’s, the Geometric Theory of Diffraction (GTD), a high-frequency 
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method, was introduced to include edge diffracted fields. The two-dimensional GTD 
model presented in [5, 6] yielded an improvement in the far-field E-plane radiation 
pattern over the approximate aperture field method. A two-dimensional model for 
the E-plane pattern, was examined by Botha et al. in [7] using an integral equation 
and the Moment Method. 

Although, the two-dimensional models perform well in the far-field E-plane pat- 
tern, they are not very accurate in predicting the H-plane pattern, especially in the 
back regions of the antenna because coupling of the diffracted fields from both the 
E- and H-plane edges is not included in the two-dimensional models. Finally, nei- 
ther the approximate aperture field method nor GTD are well suited for calculating 
aperture field distribution, VSWR, and cross-polarized patterns. 

Recent advances of the computational capabilities and the popularity of the pyra- 
midal horn antennas have encouraged the development of more accurate models with 
improved numerical efficiencies. The integral equation formulation with a Moment 
Method (MM) [8] solution has become a powerful tool in modeling complex electro- 
magnetic field problems. MM has been used to analyze an aperture on a ground 
plane [9, 10] and an aperture on a ground plane in the presence of a thin conducting 
platefl 1] . In addition, the MM has been applied to pyramidal horn antennas, both 
with and without corrugations, mounted on a ground plane[12, 13]. The presence of 
a ground plane simplifies the MM analysis. However, in most applications, the horn 
antenna is a stand-alone radiating element and is not mounted on a ground plane. 
Without the presence of a ground plane, an electric current is induced on the outer 
surface of the horn. This current has a significant impact on the radiation pattern 
at wide angles and in the back region of the antenna. 

Complete three-dimensional models have been developed for electrically small 
pyramidal horn antennas by using MM (14]. Electrically small H-plane sectoral horns 
and X-band standard-gain horns were ad so analyzed using the finite-difference time- 
domain method (FDTD) [15, 16]. High gain pyramidal horn antennas, are difficult 
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to model accurately with simple extensions of existing numerical electromagnetic 
methods, such as those described in [14, 15, 16, 17]. To model the transition from 
a relatively small feeding aperture to a much larger radiating aperture, requires the 
use of a large number of elements. Also, the electric current density on the interior 
surfaces of the horn transition becomes too complex to be modeled effectively using 
these methods. 

A full-wave stepped-waveguide model and an HFIE method to analyze both the 
interior flaring and the exterior current contributions was previously developed by 
Kuhn et a 1. [18] for conical horns. Kuhn’s formulation is simplified due to the axial 
symmetry of conical horn antennas. 

1.2 Analysis of Horn Antennas 

This part of the report presents numerical analysis methods for horn antennas with 
perfectly conducting walls. An integral equation method is first presented and then 
applied to two-dimensional models. Then a pyramided horn antenna mounted on 
an infinite ground plane is examined. This type of horn antenna mounting is often 
used on the surface of an aircraft or space-craft. Finadly, pyramidal horn antennas 
radiating in free-space is examined using a three-dimensional full wave formulation. 

The three-dimensional, full-wave formulation provides flexibility and includes all 
of the important details of a practical pyramidal horn antenna. It represents the 
first full-wave method to include the current densities on adl the conducting surfaces 
of the pyramidal horn. These are necessary to obtain the fine pattern structure in 
the regions of low-level radiation. Based on this work, a paper was submitted and 
accepted for publication in the IEEE Transactions on Antennas and Propagation 
[19]. 
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1.2.1 Two-dimensional Horn Models 
A. E-plane Pattern 

Botha et al. utilized the two-dimensional TE polarized radiation and scattering code 
by Richmond [20] to analyze the E-plane radiation pattern of a profiled horn [7]. 
The two-dimensional model can be used to simulate the principal E-plane radiation 
patterns of pyramidal horns. 

Fig. 1.1 demonstrates the E- plane geometry of a pyramidal horn. The position 
of the magnetic line source d is chosen a multiple of 4? to achieve best robustness of 
the Moment Method solution. For the two-dimensional radiation problem with TE 
polarization, the integral equation method, such as in [20], can be applied. However, 
a more simplified technique [21, 22] published by the authors of this report, presents 
a better alternative for this problem. In this simplified TE-polarized field formula- 
tion, pulse expansion and point matching were introduced. Analytical evaluation of 
diagonal elements of the impedance matrix was developed. 



L 

Fig. 1.1: E- plane geometry of a pyramidal horn antenna. 


Figs. 1.2 and 1.3 compare the normalized El-plane radiation patterns between the 
two-dimensional model and the full three-dimensional model, of 10-dB and 20-dB X- 
band standard gain horns, respectively. As shown in the figures, the prediction of the 
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E-plane radiation pattern based on the two-dimensional model provides reasonabic 
accuracy in the main beam. The agreement is better for the 20-dB standard gain 
horn because the two-dimensional model approximates a larger horn better than a 
smaller one. 
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Fig. 1.2: Normalized E-plane patterns of a 10-dB X-Band gain horn. 


B. H-plane Pattern 

The H-plane geometry of a pyramidal horn is shown in Fig. 1.4. The position of the 
electric line source d was chosen to be a multiple of to achieve best robustness 
of the Moment Method solution. This two-dimensional model was analyzed using 
integral equation formulation with TM-polarization. The simple technique outlined 
by Harrington [8] was introduced to predict the radiation pattern. 

Figs. 1.5 and 1.6 compare the normalized H-plane radiation patterns between the 
two-dimensional model and the full three-dimensional model, of 10-dB and 20-dB X- 
band standard gain horns, respectively. Predictions of the H-plane radiation pattern 






based on the two-dimensional model provide reasonable accuracy in the main beam, 
especially for the 20-dB standard gain horn. Since the coupling of diffracted fields 
from the E-plane edges are not included in this two-dimensional model, the H-plane 
patterns in the back regions were not predicted accurately. 

The advantage of the two-dimensional models is their analytical simplicity. The 
two-dimensional models predict the main beam of the radiation patterns very ef- 
fectively. The computer code for the two-dimensional models can be executed on 
a personal computer with a quick turn-around time. One of the limitations of the 
two-dimensional models is their inability to analyze the field coupling between two 
different principal planes. Therefore, back lobes are not predicted very accurately. 
Finally, the two-dimensional models can not be used to predict aperture distribution, 
VSWR, gain, and cross-polarized patterns. 

1.2.2 Pyramidal Horn Mounted on a Ground Plane 

Fig. 1.7 illustrates the geometry of a pyramidal horn antenna mounted on a ground 
plane. This problem can be separated into two parts. The first part includes the 
transition from the feeding waveguide to the radiating aperture. A complete full- 
wave approach is applied to represent this transition as a series of stepped-waveguide 
sections. This was done as shown in Fig. 1.8. Mode matching was performed by 
rigorously enforcing the boundary conditions at each step. The result is a scattering 
matrix for each step, which can be combined to obtain the scattering matrix for the 
entire transition region. 

The second part is the radiating aperture in the presence of an infinite ground 
plane. The equivalence principle was introduced to formulate the magnetic field in- 
tegral equation on the radiating aperture which combines the interior field transition 
problem with the exterior radiation. 




<j> (Degrees) 


Fig. 1.5: Normalized H-plane patterns of a 10-dB X-Band gain horn at 10 GHz. 
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Fig. 1.6: Normalized H-plane patterns of a 20-dB X-Band gain horn at 10 GHz. 
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A. Interior Horn Transitions 

Accurate analysis of waveguide transitions has been an interesting research topic in 
microwave circuit design [23, 24, 25, 26]. The available numerical approaches can 
be divided into two classes. First is a numerical solution of a system of ordinary 
differential equations [24, 27], and second is the stepped waveguide model with a 
full-wave mode-matching technique on the stepped junctions [13, 18, 23, 27]. With 
the first class, the numerical solution of the differential equations must be performed 
with a finite advancing step size. However, due to the numerical problem caused by 
evanescent modes, the taper has to be divided into several sections. The hybrid ma- 
trix of each section needs to be computed separately, then translated and combined 
into a scattering matrix. 

On the other hand, the stepped-waveguide technique uses a finite number of 
subdivisions of waveguide steps to approach the continuous horn taper. Within each 
waveguide step, the waveguide section is uniform. The scattering matrices of each 
waveguide step are related to each other by the electromagnetic boundary conditions 
on the discontinuous junctions connecting them. The combination of the scattering 
matrices of all the steps gives the total scattering matrix of the horn transition. 

As examined by [25] and [27], the numerical technique in solving the first-order 
differential equation yields results with the same accuracy as given by the stepped- 
waveguide technique, when the size of the steps is sufficiently small. The validity of 
the stepped-waveguide approximation to simulate the continuous horn transition has 
been justified in [12, 13, 18, 23, 27]. The computational effort for the two approaches 
is about the same. 

The stepped waveguide technique was employed in this project because of its 
advantage in numerical stability. Fig. 1.8 represents a typical stepped waveguide 
model of the horn transition. The continuous transition is approximated by a number 
of cascaded stepped-waveguides. For the pyramidal horn transition, each step is a 
section of rectangular waveguide. Electromagnetic field distributions in each of the 
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rectangular sections were expressed as the superposition of all possible T M t and T E z 
modes. Moreover, the TM Z and TE Z modes in each region were generated by the 2 
components of two vector potentials[22] 


F z = e Yj Mmne 30,z + B mn e }0,z )e mn (x,y) (1.1) 

mss 1 , M 
r\- Q,N 


A, = £ (C mn e- ]0 ' z - D mn e ]0 ' z )h mn (x,y) (1.2) 

m* 1 , M 
risS.N 


2 cos(/? r x) cos (f} v y) 2sin( / S r x)sin( / 3 y y) 

e mn \ x ^ y ) — _ / — 771 — 7 1 ^ranl 1 )!/) — 


0c\fab{l + Sn^ 


@ c y/ab 


A = A = y; A = -JM + % ; A ! + ft = <*> V; A- = { J 


=0 
otherwise 


Then the transverse electric and magnetic fields due to the two vector potentials 
were evaluated by the following equations [22]: 


F — 

1 dF z .1 d 2 A z 

(1.3) 

£j x — 

e dy ^ ufit dxdz 

E v = 

1 dF z .1 d 2 A z 
t dx ^ uifie dydz 

(1.4) 

H x = 

. 1 d 2 F z 1 dA z 
u\itdxdz n dy 

(1.5) 

H y = 

. 1 d 2 F z 1 dA x 

^ ui/it dydz n dx 

(1.6) 


At each stepped waveguide junction of Fig. 1.8, the boundary conditions relating 
the two slightly different-sized waveguides axe 


Afifeyllon S = y)lon s 


(1.7) 
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E?,l( x iy ) Ion A 5 

H x,l( x ,y)\on S 


0 



H%(x,y)\ 


on 


(1.8) 

(1.9) 


where 5 is the area of the smaller waveguide section, and A S is the “ring” area of 
the larger waveguide section. 

Testing the boundary conditions in (1.7)-(1.9) by the corresponding expansion 
functions (this process is equivalent to a full-domain Galerkin’s Method), uniquely 
defines the relation between the full-wave expansion coefficients on both sides of the 
junction. After some straight forward mathematical manipulations, the resulting 
transition equations for the full- wave expansion coefficients of TM Z and TE Z are 
given by 


<«£ + 32 = i: liras + Ag’) 

kl Pc Px 


(1.10) 


. n«, V wm'f - ( /W)’)l ,, , . g M) 
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y> )] u /£-(2)_p(2)\ (112) 

h rfWW* (mn mn ' 


cl! 1 - 41' = E§|f|v mn „(c<?»-i>2) 

mn jJx fJc 


(1.13) 


where 
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= ===== 
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(1.17) 


The integrals in and t/ m „« were evaluated analytically and reduced to 


n _ 4[1 + (~l) m+fc ][l + (-!)"•*•'] n ^ 
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mn “ " ^ y 4 2) 


To express ( 1 . 10)-( 1 .13) in matrix form, let 


(1.19) 


■MS) 

»(» = ( 
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J 
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k D Sl ) 
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Then, the matrix equations of (1.11)-(1.13) are 



a< 2 > + 6< 2 > 

II 

>-* 

+ 
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where Y is a matrix with corresponding elements defined by (1.10) and (1.11), and 
Z by (1.12) and (1.13), respectively. With a simple mathematical manipulation of 
the two matrix equations, the scattering matrix of the stepped- waveguide junction 
can be obtained in the form 

(/ + ZY)~\I-ZY) 2 (I + ZY)~ l Z \ 

K 2(1 + YZ)~ l Y (I + YZ)-\YZ - I) , 



a* 1 ) 


(1.24) 


The next step is to combine the scattering matrices of the stepped-waveguide 
sections. Consider two series stepped-waveguides each having the scattering matrix 
of the form 


&(’) 

a< i+l > 


)-( 


&i*» \ ( sji 

ai‘«> ) ( s<; 


S8 

c(0 

°2l 

S)( 

«w \ 

fe (.+D ) 

(1.25) 
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c(«+l) 

\ ( a (i+l) \ 


*->11 

o(*+ 1 ) 
*->21 

°12 

o(*+l) 

&22 

) ( ) 

(1.26) 


where the scattering matrices in the i-th and (i+l)-th section of stepped-waveguides 
are given by the combination of their junction scattering matrices and the corre- 
sponding propagation matrices due to the lengths of the waveguide section. These 
can be expressed as 


o(*') 

J 1X 

= [/ + Z (v) Y {v) ]~ l [l - Z^Y^} 

(1.27) 

c( 1/ ) 

°12 

= 2[I + z {l,) Y (v) Y x Z (v) D {v) 

(1.28) 

cM 

Oji 

= 2 D (u) [I + Y (v) Z^Y x Y (v) 

(1.29) 

J 2 2 

- D {v) [I + y(*')z (,,) ] _1 [y ( *' ) Z ( *' ) - I\D (v) 

(1.30) 


where the propagation matrix is a diagonal matrix expressing the phase delay 
or magnitude attenuation of all possible positively going modes in the full-wave 
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expansion of the iMh ( i-th or (i + l)-th ) section of the waveguide which has a length 
of dS u \ The element of is given by 




(1-31) 


Using (1.25) and (1.26), one can expand the two scattering matrices and obtain 
the combined scattering matrix of 


where 
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(1.34) 

(1.35) 

(1.36) 


[/ - s!; +,, s«?f 1 
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(1.37) 

(1.38) 


Finally, if all incident full-wave modes are represented by a vector a, and the 
reflected full-wave modes by a vector b, the total scattering matrix for the horn 
transition relates the two by 


c(r) a(T) 
°11 °12 


o(T) o(T) 
3-21 Ooo 


(1.39) 
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where superscripts ‘(F)’ and ‘(A)’ denote, respectively, the full-wave coefficients at 
the feed junction and radiating aperture, and the superscript ‘(F)’ denotes the total 
contribution of the horn transition. Although pyramidal horns are generally excited 
by the dominant T E\q mode from the feeding waveguide, higher order T E mn and 
T M mn modes are generated in the transition. For horn geometries with symmetrical 
flaring in both dimensions, only those modes which have m = 1,3,5, 7,..., M and 
n = 0,2,4, 6 ,..., N (n = 0 for TE modes only) are generated in the transition. The 
coupling between T E mn and T M mn modes is automatically included in the full-wave 
analysis of the stepped junctions. 

There are two factors which affect the accuracy of the stepped-waveguide ap- 
proximation. First is the size of the steps and second is the number of terms in the 
series ( M and N in ( 1 . 1 )-( 1 .2)) . Past experience demonstrated that convergence is 
achieved by limiting the maximum size of the steps to A/32. The choice of M and N 
should be based on the flare Jingle and length of steps. For horn antennas mounted 
on a ground plane, the choice of M and N is less critical. A choice of M = 5 and 
N = 4 is sufficient to predict accurately both the antenna patterns and the aperture 
field distributions. 

B. Integral Equations 

The radiation from the aperture into the half-space can be analyzed using an integral 
equation formulation and a Moment Method solution [8, 22]. Since the aperture is 
mounted on an infinite ground plane, only a magnetic field integral equation equation 
is needed to solve the aperture discontinuity. This equation is 

2H“‘(M) = on the aperture (1-40) 

where field components with ‘ext’ denote the fields on the half-space side of the aper- 
ture, and ‘int’ the fields on the horn side of the aperture. The equivalent magnetic 
current density is given by 
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M(a< A > + b ( ^) = -n x E in ‘(a (>l) , b^>) (1.41) 

where n is the unit vector normal to the aperture, and E ,nt (a^, b*' 4 ') is the electric 
field internal to the aperture. The Moment Method solution of the integral equation 
is a simplified procedure from the field problem of a pyramidal horn antenna in free 
space presented in the next section. 

Recently, a technique was introduced by Williams et al. [28] to simulate ex- 
perimentally radiation of an antenna on an infinite ground plane using available 
measurements from the same antenna mounted on a finite-size ground plane. The 
foundation of the technique presented by Williams et al. is based on the fact that 
the scattered electromagnetic fields on the front and back sides of an electrically thin 
ground plane are symmetric. The measured radiation pattern in the back-side of the 
finite-size ground plane is used to correct the front-side of the radiation patterns by 
subtracting off the diffracted fields. Therefore, the corrected front-side radiation pat- 
terns represent the radiation patterns of the antenna mounted on an infinite ground 
plane. 

This “subtraction” technique was applied in this project to convert measured 
radiation patterns of a pyramidal horn mounted on a finite-size ground plane to 
corresponding radiation patterns of a pyramidal horn mounted on an infinite ground 
plane. The radiation patterns for the latter case was computed numerically and 
compared with the corrected experimental measurements. For our experiment, a 
20-dB standard gain horn antenna (A = 4.87", B = 3.62", L = 10.06", a = 0.9" and 
b = 0.4") was mounted on a 3 x 3 feet planar aluminum sheet ground plane. 

Figs. 1.9 and 1.10 compare the measured and computed E- and H-plane radiation 
patterns of the 20-dB standard gain horn mounted on an infinite ground plane. The 
measured padiation patterns from a 3 x 3 finite-size ground plane were corrected in 
Figs. 1.9 and 1.10 to correspond to the infinite ground plane using the “subtraction” 
technique in [28], As illustrated by both figures there is very good agreement between 
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the experimental and computed radiation patterns. 

The results obtained from this study lead us to conclude that the analysis of 
pyramidal horn antennas mounted on a ground plane can be performed accurately 
and efficiently using the stepped-waveguide and Moment Method procedure. 



Fig. 1.9: El-plane radiation patterns of a 20-dB X-Band standard gain horn at 10 GHz, 
mounted on an infinite ground plane (A = 4.87", B = 3.62", L = 10.06", a = 0.9" and 
b = 0.4"). 
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Fig. 1.10: H- plane radiation patterns of a 20-dB X-Band standard gain horn mounted on 
am infinite ground plane. 
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1.2.3 Pyramidal Horn Antennas in Free Space 

Fig. 1 . 1 1 illustrates a pyramidal horn antenna in free space. The field problem in 
this case also can be separated into two parts. The first part includes the transition 
from the feeding waveguide to the radiating aperture. The second part includes the 
aperture radiating into the free space. 

The interior transition from the feeding waveguide to the aperture can be ana- 
lyzed using the stepped-waveguide model as outlined previously for the analysis of 
pyramidal horns mounted on a ground plane. However, in order to obtain accurate 
sidelobe levels of the antenna pattern, a sufficient number of terms in the series must 
be included because higher order modes at the aperture contribute significantly to 
the radiation pattern in the back regions. Significantly more terms axe needed for this 
analysis than that in [12, 13] for pyramidal horns mounted on a ground plane. An 
empirically derived formula for determining M and N is the nearest higher integer 
of 

(M, N) = — + 1.5 (1.42) 

where A and B are the dimensions of the horn aperture. When analyzing a high- 
gain horn antenna, the required number of modes determined by (1.42) is large, 
and if used throughout the matching process, the computation could become very 
inefficient. To reduce the computation time, a variable number of modes was used 
along the transition. At each step, the number of modes was determined using the 
same criterion given for the aperture in (1.42) except A and B are replaced by the 
step dimensions. Therefore, only a few modes are required near the feed, whereas 
the number of modes approach (1.42) as the computation proceeds toward the horn 
aperture. This process preserves accuracy, and for the geometries considered in this 
research reduces the computation time by eight compared to that of a constant 
number of modes. 

Due to the absence of the ground plane on the aperture, the modeling of the 
antenna aperture is more complicated. The following sections outline a hybrid field 
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Fig. 1.11: Geometry of the pyramidal horn antenna in free-space. 



22 


integral equation formulation for the pyramidal horn aperture radiating into free 
space. 

A. The Hybrid Field Integral Equation 

Fig. 1.12 presents the general problem of a radiating aperture on a conducting body. 
The fields internal to the aperture are represented by the full-wave vectors a (the 
incident mode) and b (the reflected modes). To relate the internal and external fields 
on the aperture, Love’s field equivalence principle [22] was introduced. The radiating 
aperture was replaced by a sheet of perfect electric conductor with a magnetic current 
density given by 

M(a ( ^ + b (>l) ) = -n x E ,nt (a (>1) ,b (> ‘ ) ) (1.43) 

where n is the unit vector normal to the aperture, and E ,nt (a^), b^) is the electric 
field internal to the aperture. Equation (1.43) insures the continuity of tangential 
electric fields across the aperture. The magnetic current density M is radiating in 
the presence of the closed conducting surfaces of the horn. An electric surface current 
is then induced to maintain zero tangential electric field on the external surfaces of 
the horn. The continuity of the tangential magnetic field across the aperture and 
the boundary conditions on the external surfaces require that 

H“*(M) + HJ^(J) = HJ”J*(a (jt) ,b (>l) ) on the aperture (1.44) 

E^(J) + E^*(M) = 0 on the external surfaces and the aperture (1.45) 

where field components with ‘ext’ denote the fields on the free-space side of the 
aperture. The boundary conditions of (1.43)-(1.45), derived from Love’s field equiva- 
lence principle , satisfy the uniqueness theorem and define an equivalent to the actual 
electromagnetic problem. 
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B. Moment Method Solutions 

In the Moment Method solution of the hybrid field integral equations of (1.43) and 
(1.45), subsectional quadrilateral roof-top patch modes were chosen as both expan- 
sion and testing functions for the electric current distribution on the exterior surfaces 
of the horn. To exploit the Toeplitz property of the impedance matrix elements, the 
magnetic current on the aperture was also expanded and tested with the same roof- 
top patch modes as those of electric current modes on the aperture. Thus, the 
impedance and admittance matrix elements for the integral equation solution on the 
aperture were related and were not computed twice. A more detailed explanation of 
such this process can be found in [28]. 

Since the aperture magnetic current density is related to the full-wave modes 
and b^), a conversion matrix is introduced to transform the roof-top patch modes 
into eigen modes of the aperture field in (1.43) 

Nj N m N a 

J = M = £M,m t = X>,£K J Pf (1-46) 

t=i i=i » ;=i 

where P/ and are the subsectional roof-top patch modes for the i-th electric 
and j - th magnetic current density modes, respectively; Nm is the total number of 
full- wave expansion modes of both TE and TM in the aperture. Na is the number 
of roof-top patch modes representing the aperture magnetic current. Nj is the total 
number of roof-top patch modes representing the electric current on the entire outer 
surface of the horn including the aperture. [V^j] is the conversion matrix from the 
distribution functions of the aperture modes m,- expressed by e mn (x, y) and h mn (x, y) 
in (1.1) and (1.2) to the roof-top patch modes. Testing (1.45) with P/ leads to 

-<P/,E“(J)>=<P/,E“(M)> (1.47) 


Substituting the expansion representation of J and M into (1.47), yields a matrix 
with elements given by 

Z* - -i < Pf.ESUP/) > = -1 X L Pf ' E£(P/)<M»j ( 1 . 48 ) 

T 1 7 JSi JSj 
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g H =<P^E'"(P")> = / f Pi E^f)dHds, (1.49) 

Js k Js, 

A normalization factor of I /77 is introduced in computing the matrix elements of Z tJ 
in (1.48) to reduce rounding errors in combining large matrices, which is very helpful 
in maintaining the numerical stability of the MFIE in (1.44). Efficient and accurate 
algorithms are developed for evaluating the four-fold integral of (1.48) and (1.49). 
Descriptions of the algorithm can be found in Appendix A. 

Since the expansion and the testing modes are the same, the impedance matrix 
Z is a complex symmetric matrix. It can be filled and decomposed with only half of 
the CPU time needed for the full-matrix system [29, 30]. The solution of (147) is 
given by 

J = -Z -1 UM (1.50) 

V 

U = QV‘ (1.51) 

where the superscript ‘t’ signifies transpose. J is the vector representation of the 
subsectional roof-top patch modes, and M represents the aperture field modes. The 
transformation from Q to U in (1.51) reduces the number of right-hand side solution 
from Na to Nm ( Nm is generally a small fraction of Na ). Again, testing (1.44) with 
m; , the MFIE can be expressed by 

< m,-,H“‘(M) > + < m,,H“‘(J) >=< m,-,H£,(a< A >,b< A >) > (1.52) 

Using duality for the reaction of aperture magnetic current modes[28], reciprocity 
for the reactions between the aperture magnetic mode and the outer surface electric 
modes, and the relation of (1.48), (1.52) can be written as 

b (i4) = S {A)n sS A) = [Y (a) + Y‘ n *] - 1 [Y‘ n< - Y (i 4 ) ]a ( ^ (1.53) 
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where I is the identity matrix, Y^, and Y* mt ) the aperture admittance matrices 
defined by 

YM)( a (A) +b (A) ) = [VZ (>l) V t + U‘Z- 1 U][<m 1 ,M(a M) + b (A) ) >] (1.54) 

Y in ‘(a (>l) - b {A) ) = [< - b (A) ) >] (1.55) 

where Z (A) is a subset of Z for electric current modes on the aperture. Using (1.53) 
and (1.39), the incident aperture field distribution coefficients (denoted by a vector 
a^*) and the reflected wave at the feed (denoted by a scalar coefficient 6^) are 
given by 

SL tA > = (I-S'rsK’J-'S'TM? (1.56) 


iff = [SST* + Sg ) S< J ‘*(I - Sg 1 S , ' ,) )- 1 S^ l ]o ( i5 1 (157) 


where is the magnitude of the incident TE\q mode in the feeding waveguide. 
Equations (1.57), (1.56), (1.53), and (1.50) specify the antenna reflection coefficient 
(or return loss), the aperture magnetic current distributions, and the electric current 
distributions on the exterior surfaces. The radiated electric field can be determined 
from the current expansions and the gain computed from 


G(9,4) — 101ogi O 


|E(J,M) + E(M,M)r 


(1.58) 


where f c is the cut-off frequency of the feeding waveguide, and / is the operating 
frequency. 


1.3 Results and Discussions 

Five X-band high-gain pyramidal horns were chosen for this analysis. Table 1.1 
gives the dimensions and typical computed data for a maximum segment length of 
0.15A for 10- and 15-dB standard gain horn, and 0.2A for the other three horns in 
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the quadrilateral roof-top expansion of the exterior surfaces at 10 GHz. The CPU 
times are typical for an IBM R6000-350 workstation. 

Table 1.2 lists the gains and VSWR’s obtained by the approximate method and 
the technique of this report. As expected, differences are more obvious for 10-dB gain 
horn since the approximate method [4] does not work very well for such a small horn. 
Table 1.3 compares approximated, measured, and predicted gains and VSWR’s for 
the 20-dB X-band standard gain horn antennas at three different frequencies. Note 
that the gains listed in the tables have almost a constant 0.2 dB difference between 
the moment method solution and the measured data. The same amount of difference 
is observed between computed and measured gains for the other two large X-band 
square aperture horns. It should be mentioned that the same 20-dB standard-gain 
horn is used as the calibration antenna. Therefore, if there is any inaccuracy in the 
calibration data, all measured data sets are affected. The agreement between the 
computed and measured VSWR’s is excellent. 

Table 1.1: Typical data of pyramidal horn antennas analyzed 



Pyramidal 

A 

horn 

B 

dimensions 

L 

EFIE Matrix 
size(# of rows) 

CPU 
(in hrs) 

10-dB Standard-gain horn 

1.58" 

1.15" 


820 

0.38 

15-dB Standard-gain horn 

2.66" 

1.95" 

5.46" 

1600 

0.61 

20-dB Standard-gain horn 

4.87" 

3.62" 


4300 

1.8 

5- inch square horn 



1 

5700 

2.6 

7-inch square horn 



12.1" 

9600 

11.3 


Figs. 1.13 to 1.14 compare the E- and H-plane patterns obtained from this 
method and the approximate method for the two smaller X-band horns, respectively. 
The approximate patterns are computed using the method outlined in Chapter 12 
in [4], except that the free-space wave impedance in (12. Id) is replaced by the guide 
wave impedance at the aperture. As expected, the agreement in Fig. 1.13 is not as 
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Table 1.2: Comparison of VSWR’s and gains of 10- and 15-dB standard gain horns 




8.2 

GHz 

10.3 

GHz 

12.4 

GHz 



VSWR 

Gain 

VSWR 

Gain 

VSWR 

Gain 

10-dB 

Approximate 

N/A 

8.98 dB 

N/A 

11.13 dB 

N/A 

12.81 dB 

Horn 

This Method 

1.18 

9.75 dB 

1.17 

11.63 dB 

1.20 

13.48 dB 

15-dB 

Approximate 

N/A 

13.83 dB 

N/A 

15.83 dB 

N/A 

17.42 dB 

Horn 

This Method 

1.11 

14.23 dB 

1.14 

15.94 dB 

1.10 

17.58 dB 


Table 1.3: Comparison of VSWR’s and gains of the 20-dB standard gain horn 



9 GHz 
VSWR Gain 

10 GHz 
VSWR Gain 

11 GHz 
VSWR Gain 

Approximate 
This Method 
Measured 

N/A 19.77dB 
1.082 19.98dB 

1.10 19.72dB 

N/A 20.59dB 
1.057 20.63dB 
1.06 20.46dB 

N/A 21.31dB 
1.031 21.46dB 

1.04 21.24dB 


good as in Fig. 1.14. 

Fig. 1.15 compares the computed, approximated, and measured E- and H-plane 
patterns for the 20-dB standard gain horn. Since the horn is much larger than the 
previous two horns, the approximate method compares well with the measurements in 
the first few lobes of the patterns. However, the patterns predicted using the moment 
method code compare much better with the measured patterns. The improvement 
over the approximate method is more evident in the back region of the antenna. 

In 1965, Russo et ai. [5] presented an investigation of using different wadi thick- 
nesses to model the horn and their effect on the E-plane patterns of pyramidal horns 
using GTD. Fig. 1.16 presents our comparisons of the E- and H-plane patterns for 
the 5-inch square X-band horn at 10.0 GHz with and without modeling the aperture 
wall thickness. This comparison agrees with Russo’s conclusion for the E- plane ra- 
diation patterns. However, for the H-plane patterns, the predicted patterns of the 
thin wall model do not agTee in the back region with the measured patterns as well 
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as those of the thick wall model. The thick wall model accurately predicts the fine 
ripple structure in the back region while the thin wall model does not. The aperture 
wall thickness affects the fine ripple structures of the H-plane patterns. The addition 
of the outer surface of the feeding structure in the EFIE has only a negligible effect 
on the patterns of the horn antennas. 

Another advantage of the solution implemented in this research is the ability 
to examine the aperture field distribution. Figs. 1.17 to 1.19 represent computed 
aperture field ( E y and H z components) distributions of the three X-band pyramidal 
horns at 10 GHz, respectively. Contrary to the approximation that the aperture 
fields are basically T Ei 0 with parabolic phase fronts, the distributions presented 
here demonstrate a more complicated shape. The deviations from the approximate 
aperture distributions are more evident for H x . Figs. 1.17 to 1.19 also illustrate 
that as the electrical size of the aperture become larger, the amplitude distributions 
become more complicated, but the phase distributions approach a parabolic phase 
front as assumed by the approximate method. 

As presented in the comparisons between computed and measured results, the 
full-wave HFIE method has demonstrated an excellent accuracy both for small find 
large pyramidal horn antennas. One of the requirements of this method in solving 
large pyramidal horns is computer memory. In some workstations such as the IBM 
RS6000/350, the fast data transferring capability between the memory and the hard- 
disk enables the out-of-core memory matrix solver for a symmetric complex system 
to run nearly as fast as the in-core solver. If the symmetric property of the electric 
current distribution on the exterior surfaces is utilized, the memory can be reduced 
to 1/8 while the CPU time cam be decreased to about 1/4 of the values listed in 


Table 1.1. 
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Fig. 1.14: Comparison of E- and H-plane patterns for 15-dB standard gain horn at 10 GHz. 
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Fig. 1.15: Comparison of E- and H-plane patterns for 20-dB standard-gain horn at 10 GHz. 
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Fig. 1.17: Aperture fields of 10-dB X-band standard-gain horn at 10 GHz. 


Aperture fields of 10-dB standard gain horn at 10 ClHa. 

(a), (b) Magnitude and phase distributions of 

(c). (d) — Magnitude and phase distributions of H_ 









36 



Fig. 1.19: Aperture fields of 7-inch horn at 10 GHz 


Aperture fields of 7 — inch square horn at 10 GHz. 

(a), (b) Magnitude and phase distributions of 

(c). (d) Magnitude and phase distributions of H 
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1.4 Conclusions 

A full-wave, hybrid field integral equation method has been developed to analyze 
electromagnetic horns of pyramidal configuration. The model includes the current 
induced on the exterior surfaces and has been verified by comparison of computed 
and measured data. For accurate results, the following guidelines are suggested: 

• include a sufficient number of higher order modes in the full-wave analysis. The 
required number of modes is determined using the empirical formula (1.42). 

• limit the size of the stepped discontinuity to less than A/32 when approximating 
continuous horn transition. 

• limit the largest segment size of the roof-top patch to less than 0.2A to ensure 
an effective and accurate solution of the electric field integral equation on the 
outer surfaces of the pyramidal horn antenna. 

• include the wall thickness in the segmentation of the outer surface. 

The research in this project demonstrated that the approximate method found 
in most antenna books is not very accurate in predicting gains and patterns of the 
pyramidal horns of small electrical size. However, as the electrical size gets larger, 
the approximate method becomes more accurate in predicting gains and the first few 
sidelobes of the fax-field radiation patterns of a pyramidal horn. The aperture field 
distributions of pyramidal horns are much more complicated than those predicted by 
the approximate formula. However, the phase distribution approaches a parabolic 
distribution for apertures of large electrical sizes. 

Although a straight forward application of this method is computationally inten- 
sive, a widespread application of high-power workstations, such as the IBM RISC6000 
platforms, is making it more effective to serve as a computer-aided analysis and design 
tool. Although the present development has concentrated on pyramidal horn anten- 
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nas, the model can be extended to analyze wide-band dual-ridged, and quadruple- 
ridged horns. 


CHAPTER 2 


MOMENT METHOD ANALYSIS OF HORN ANTENNAS WITH LOSSY WALLS 
2.1 Introduction 

Most horn antennas with simple designs, including the pyramidal horn antenna, 
have asymmetric E- and H-plane radiation patterns and relatively high cross-polar- 
ization levels. These two properties become undesirable when the horns are used 
as feeding devices for reflector antennas. Many complex designs were developed to 
overcome these two undesirable features[l]. One was to introduce corrugations in the 
transition from the feeding waveguide to the radiation aperture. The corrugations 
on the interior surfaces of horns successfully suppress the sidelobes of the E- plane 
radiation pattern, maintain symmetric E- and H-plane radiation patterns, and reduce 
the cross-polarization [31, 32]. However, due to the high cost and heavy weight of 
the corrugations, application of corrugated horn antennas is still limited. The search 
for the substitution of the corrugations has become an interesting topic for many 
years [33, 34]. 

A waveguide with lossy wall was undesirable until Knop et a i. discovered its 
antenna beam shaping capability in the early 1980’s [33]. Independently, Lee et 
aJ. also noted the same field tapering in a circular waveguide by introducing a 
thin coating of lossy material on the waveguide[35, 36]. Knop et a 1. also applied 
this finding to shape the radiation pattern of a large conical horn antenna [37, 38]. 
Wang et a 1. introduced lossy magnetic coating to produce the same beam shaping 
results on pyramidal horns [39, 40]. Ghobrial et a 1. introduced lossy walls to suppress 
the cross polarization of a pyramidal horn used as a feeding antenna for a paraboloidal 
reflector [41]. These early experimental studies have demonstrated the effectiveness 
of lossy material coatings in controlling the radiation patterns of horn antennas. 
Although a few simplified methods were presented to provide guidelines in using 
lossy material coatings on the walls of horn antennas [36, 37, 38, 39, 40], a rigorous 
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analysis for horn antennas with impedance walls has not been reported. 

The full-wave, hybrid integral-equation method presented in this report provides 
an analytic solution to the pyramidal horn antennas with or without the imped- 
ance walls. The validation and efficiency of this analytical method is important in 
the application of lossy material coatings to control the radiation pattern of horn 
antennas. 

2.2 Eigen Expansions 

The coating of lossy materials on the walls of the horn antenna introduces a change 
of the electromagnetic field distribution along the horn transition. It is important to 
note that maintaining a symmetric coating is very important because asymmetry in 
the coating will introduce significant cross polarization, as noted in [41]. 

Fig. 2.1 shows the geometry of an E-plane coated horn transition and its stepped- 
waveguide model. One difference between the coated and uncoated stepped waveg- 
uide sections is the eigenfunctions in the full-wave expansions. To find the new 
eigenfunctions, the new eigenvalues for each individual eigenfunction based on the 
boundary conditions across the surfaces of the coating material must be found. Such 
process can be analyzed by using the Transverse Resonance Method (TRM) [22] 
when coatings are applied on the top and the bottom walls, and TM y and TE y 
modes jure used to represent the field distribution. By using the TRM technique, the 
cross-sectional wave constants can be obtained by solving the following two complex 
transcendental equations for TE y and TM V modes respectively, are obtained: 

— tan (P e h) + Pm tan(A*d) =0 p 2 - Pm = w 2 /i 0 c 0 (/i r e r - 1) (2.1) 

(■d 

— COt(P m h) + Pm COt (Pond) = 0 P m 2 - P^ 2 = U) 2 fi 0 e 0 (fi r t T - 1 ) (2.2) 

Ed 

where (&,/?,*.), and ( P m ,Pom ) are the complex propagation constants inside and 
outside the coating material for T E y and T M y modes, respectively. 
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The solution of (2.1) and (2.2) requires a complex root solver to obtain all the 
roots in the complete eigenvalue sequence. Such a complex root solver has been de- 
veloped for this project. Once the complex eigenvalues are provided, the scattering 
matrix for each individual coated waveguide step can be found following a similar 
procedure as the one applied to the perfectly conducting waveguide step. The tech- 
nique of combining the scattering matrices of the coated- waveguide steps is identical 
to that of the uncoated case. 

When the modified electromagnetic field distribution (due to the presence of 
the lossy coating) is computed using the complex eigenfunctions, the results are 
analytical. The effect of impedance surfaces is therefore computed more accurately 
using the complex eigenfunctions than by using an approximate impedance boundary 
condition. 

2.3 Low-Loss Design Using Impedance Surfaces 

Investigation of using impedance surfaces in horn antennas for radiation pattern 
control have been conducted by many researchers [36, 38, 40]. The introduction of 
lossy material on the inner walls of the horn antenna presents a good alternative 
to corrugations. However, the coating of lossy material introduces ohmic losses 
ranging in gain losses from 0.5 dB in [38], 10 dB in [36], and 0.8-2. 7 dB in [40]. Our 
objective in applying impedance surfaces to horn antennas is to control the horn 
antenna radiation pattern with minimum ohmic loss. This can be achieved through 
a comprehensive study of the effect of the coating on the the radiation patterns of 
horn antennas. 

Based on the stepped-waveguide mode outlined in Fig. 2.1, the ohmic loss is 
mainly contributed by the imaginary part of the propagation constant of the dom- 
inant mode (HEu) along the horn transition. Further studies of the eigenvalue of 
the dominant mode (HEu) suggest that the ohmic losses depend upon the material 
properties, the thickness of the coating, and the size of waveguide cross section [42]. 
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The following design guidelines hightlight the importance of each factor: 

• As was noted in [36], the ohmic losses decrease monotonically as the cross- 
sectional size increases because the electromagnetic fields are expelled from the 
lossy material when the cross-sectional size of the horn transition is large. This 
explains why the antenna in [38] has significantly smaller ohmic losses than 
others. One can reduce the ohmic losses by not applying the lossy coating on 
the small cross-sectional portion of the horn transition while applying a robust 
coating on the large cross-sectional portion of the horn. 

• For a large cross-sectional size antenna, the choice of the material is less critical. 
For a smaller size antenna, choice of the lossy material is more important to 
achieve the needed HE\\ mode in the horn transition. A good material for 
such an application is one having a high magnetic loss. 

• There exists an optimum thickness of the lossy material coating for low-loss 
implementation of the HE\\ mode. 

Both analytical and experimental results presented in a later section demonstrate 
a significant reduction of ohmic loss by following the above design guidelines. 

2.4 Horn Transition with Impedance Walls 

In previous reporting periods, we have presented the stepped-waveguide model for 
analyzing the lossy impedance coating on the interior surfaces of the horn transition. 
Figure 2.1 represents the stepped waveguide model for the analysis of a partially 
coated horn transition. Unlike the perfectly conducting part of the horn transi- 
tion, the stepped-waveguide model for the lossy material coated section of the horn 
transition exhibits a much more complicated boundary conditions at the stepped 
discontinuity. 

Fig. 2.2 shows three types of interfaces encountered in the stepped model. The 
difficulty involved in an accurate analysis of the lossy coating can be discussed in 
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two parts. First, to avoid a high gain loss, the thickness of the lossy material coating 
should be relatively thin compared to the cross-sectional sizes of the horn transition. 
Therefore, to effectively control field tapering, the coated lossy material should have 
a relatively high dielectric constant and/or a high permeability[36, 39, 40]. The 
analysis of the field tapering can be carried out by using eigen modes provided by 
the partially-filled waveguide techniques [22]. An accurate expansion of the field 
distributions controlled, by a small part of the discontinuity at the edges of the cross 
section, requires many eigen modes. This requirement presents a challenge in the 
numerical implementation. Second, the stepped-waveguide model artificially creates 
media discontinuities at the stepped junction which could act as a false control on 
the field distributions. 

To investigate these problem, three different techniques in modeling the coating 
discontinuity were studied in this project: 

• Solving the stepped discontinuity exactly as shown in Figure 2.2. This tech- 
nique is an exact boundary condition solution, if a real stepped discontinuity 
at the interface actually existed. However, due to the incompleteness of eigen 
modes and the artificially created material discontinuities, this method does 
not work very well. The predicted El-plane pattern presented in the previous 
report showed some discrepancies even in the main lobe of the pattern. 

• Approximating the material interface by assuming a continuous material tran- 
sition (there is no material-air interface as shown in Fig. 2.2 for the stepped 
discontinuity). By this assumption, the eigen values of eigen modes on both 
sides of the stepped discontinuity are assumed to have their original values. 
This method average out some of the artificial discontinuities on the interface. 
As a result, the computed patterns have demonstrated a better agreement with 
measured data. 

• Hybrid Fourier transform method. This method applies a Fourier transform to 
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Fig. 2.2: Stepped discontinuity of in a stepped waveguide model. 
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expand the field distributions of eigen modes on both sides of the stepped dis- 
continuity. By testing the transformed field distributions with different weight- 
ing functions at the discontinuity, the full-wave coefficients are related. Once 
the discontinuity is resolved, the field distributions are transformed back to the 
eigen mode expansion in the waveguide section to account for the propagation 
factor. This method provides the best agreement with measurements among 
the three methods considered. We do not yet have a sound explanation for 
this result. However, it looks like this method averages out the error more 
effectively than the other two. 

2.5 Results and Discussions 

Based on the hybrid Fourier transform stepped discontinuity algorithm, numerical 
simulations were carried out to predict the El-plane radiation patterns for the 20-dB 
X-band standard gain horn with 2cm, 5cm, and 10cm Northrop Nitrile lossy magnetic 
material coating from the radiating aperture. The Nitrile material has a measured 
relative permittivity of 18.75 — jO.l and a relative permeability of 1.55 — j 1.85 at 10 
GHz. The thickness of the material is 44 mils. However, considering the tilting effect 
when the material is coated, the effective thickness is taken to be 44.41 mils. In Figs. 
2.3 through 2.5, the predicted patterns of the horn without a material coating (PEC 
wall surfaces) are included to demonstrate the effect of lossy material coatings on 
the control of antenna radiation patterns. 

Fig. 2.3 compares the measured and predicted El-plane radiation patterns for the 
20-dB X-Band standard gain horn with 2cm of lossy magnetic coating. We observe 
a good agreement in the main lobe and small discrepancies in the back region. This 
demonstrates that the Fourier transform technique can be used to analyze the lossy 
coatings with reasonably good accuracy. Fig. 2.4 represents the El-plane patterns 
for the same horn with 5cm of lossy magnetic material coating. It exhibits the same 
agreement as that of Fig. 2.3 except that the measured gain is a about 0.7 dB lower 
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than predicted. Fig. 2.5 compares the E-plane patterns for the horn with 10cm 
lossy magnetic coating. In this figure, more discrepancies between the predicted and 
measured patterns are observed. It is believed that the discrepancies are caused by 
the following two sources: 

• The inaccuracy of the stepped waveguide model in the lossy horn transition. 
A thin lossy material coating on the E-plane surfaces of the pyramidal horn 
introduces an abrupt discontinuity for the electric field of the dominant mode. 
This discontinuity is more severe for the Nitrile material which has a very high 
dielectric constant. Since the material is very thin compared to the horn transi- 
tion, it is very difficult for any numerical method to model very accurately and 
efficiently. Such a disagreement becomes more evident when a larger portion 
of the horn transition is coated by the lossy material. 

• The imhomogeneous lossy material. The lossy materials used in the experimen- 
tal verification may not be homogeneous. Therefore, the material constants 
used in the analysis, which is measured from samples of the material, may not 
have accurate representation of the material actually used in coating the horn 
surfaces. 

Figs. 2.3 to 2.5 also demonstrate that the application of the lossy materials on 
the E-plane walls has a strong control on the E-plane radiation patterns. The length, 
the thickness, and the properties of the material coating control the sidelobe level of 
the E-plane in a very sensitive manner. Therefore, there exists an optimum design 
of material coating for beam shaping without a significant loss of the antenna gain. 
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Fig. 2.3: Comparison of E-plane patterns for standard gain horn at 10 GHz with 2cm of 
lossy Nitrile material coating. 
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2.6 Conclusions 

A full-wave, hybrid field integral equation method has been developed to analyze 
pyramidal horn antennas with lossy walls. The method can be used to analyze 
parti ally- coated pyramidal horn antennas. The presented technique achieves a rea- 
sonably good agreement between theoretical predictions and experimental data. This 
analytical technique represents the first full-wave technique for the analysis and de- 
sign of horn antennas with lossy walls. Based on the the eigen mode analysis, guide- 
lines are provided to design pyramidal horn antenna with lossy walls to shape the 
antenna pattern with minimal loss of the antenna gain. 




CHAPTER 3 


CONTOUR PATH FDTD ANALYSIS OF HORN ANTENNAS 
3.1 Introduction 

In addition to the full-wave Moment Method analysis, an alternative approach based 
on the Finite-Difference Time- Domain (FDTD) method was also applied in the anal- 
ysis of pyramidal horn antennas with and without composite E-plane walls. Because 
of the flared surface of the pyramidal horns in the E- and H-planes, some augmenta- 
tions of the basic FDTD method had to be made to model accurately and efficiently 
the surface of the horns. 

Pyramidal horns with and without lossy walls were modeled by applying a contour 
path FDTD method to model the antenna flared surfaces. The idea behind the 
contour path approach is to distort the FDTD grid where it is only needed (near the 
antenna flared surface) and use the basic FDTD algorithm everywhere else in the 
numerical domain. Modified finite-difference equations are applied over the distorted 
grid. The equations are obtained by applying Maxwell’s integral equations over the 
distorted part of the grid and making assumptions regarding the distribution of the 
electric and magnetic fields in the distorted grid. 

In this report some familiarity with the basic FDTD method is assumed. The con- 
tour path FDTD method is then introduced based on the integral form of Maxwell’s 
equations. This method is then used for the modeling of metallic pyramidal horn 
antennas. In the last part of the report, the contour path FDTD method is applied 
to the modeling of pyramidal horn antennas with lossy inner El-plane walls. In all 
cases, numerical results jure compared with available measurements demonstrating 
the accuracy of the contour path FDTD approach. 
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3.2 FDTD Method Based on the Integral Form of Maxwell’s Equations 
Maxwell’s equations can be written in integral form as 

jfH-dl = j t j J'tE ■ ds + J J'<rE ■ ds (3.1) 

jfEdl = -2-J J^H-ds- j f' P Hds (3.2) 

where (3.1) and (3.2) are referred to as Maxwell’s integral equations as derived from 
Ampere’s and Faraday’s laws, respectively. E is the electric field in (V/m), H the 
magnetic field in (A/m), e the electrical permittivity in (F/m), a the electrical con- 
ductivity in (S/m), p the magnetic permeability in (H/m), and p the magnetic re- 
sistivity in (D/m). The electrical permittivity t and the magnetic permeability p of 
composite materials are often provided in the form e = d — je" and p = p' — jp", re- 
spectively. The electrical permittivity and magnetic resistivity, a and p, respectively, 
are then represented by: 

<7 — U)t" — ( 3 - 3 ) 

p - up" = upop", (3.4) 

where cq and po are the electrical permittivity and magnetic permeability of free 
space, respectively, and u = 2ir /, where f is the frequency. Using the above two 
expressions the imaginary parts of electrical permittivity and magnetic permeability 
are lumped into <r and p, respectively. 

Equations (3.1) and (3.2) are implemented in the conventional FDTD method 
over spatially orthogonal contours referred to as Ampere’s and Faraday’s contours, 
respectively [43]. Fig. 3.1 illustrates an example of Faraday’s and Ampere’s con- 
tours. Application of (3.1) and (3.2) over a cubic grid proceeds by assuming that the 
dimensions of the cell are small enough so that the electric and magnetic fields can 
be considered linear functions of position, and their values at the midpoints of each 
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Fig. 3.1: Ampere’s and Faraday’s contours for implementing Maxwell’s equations in in- 
tegral form. 
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side of the contour equal to the average value of the field along that side. Under 
these assumptions, the resulting finite-difference equations for advancing the electric 
and magnetic fields are exactly the same as those obtained by applying Yee’s algo- 
rithm [43]. This formulation is second-order accurate because the field components 
are assumed linear over the cells. 

Consider the application of (3.2) over the Faraday’s contour shown in Fig. 3.1. 
This results in 

jf E dl = ~jJ l pH x dydz (3.5) 

where sx is the area enclosed by contour cj. It was assumed that ci encloses a 
free-space cell for which p = 0. The positive direction for ci is taken in the counter- 
clockwise direction in accordance with the right-hand rule. Evaluation of (3.5) pro- 
ceeds by assuming that the dimensions of the cell are small enough so that the fields 
can be assumed linear functions of position. Consider first the left-hand side of (3.5) 
which involves four separate integrals, one integral per each side of contour cj. As- 
sume that the electric field distribution on the left side of contour Ci to be of the 
form 

E, = E,(i,j,k + ^)+A-z (3.6) 

where z is measured with respect to the center of the cell and A represents the slope 
of the electric field distribution. E t (i, j, k + 1) represents the value of the field at the 
center of the left side of contour Cj. The lower left corner of the cell is denoted as 
cell (j, j, k). Thus the contribution from the left side of Ci is 

/ L E * dz ~ ~ / [ E * (*’■*’ k + + A ' Z ] dz = ~ Ex (*’■*’ k+ \)' Az ^ 3 ' 7 ) 

where the negative sign was used because the electric field component points in the 
opposite direction (clockwise) from that of contour C\. Similar integrations can be 
carried out over the rest of the contour to obtain: 

£E-dl = [f?,(i,j + l,Jfe + i) -E,(i,j, * + 5 )] A* 
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+ E v + 2 ’^) — E y (i,j + 2’^+ l)] ^ y 


The x component of the magnetic field must now be integrated over the area enclosed 
by contour C\. Assume that H x is linear within Si and its distribution over the cell 
can be expressed as: 

1\ „ 


H x = H t (i,j + ±,k + ^+B-z + C-y 


where y and z are measured with respect to the center of the cell. H x (hj + ^ + j) 

represents the value of the field at the center of the cell and the constants B and 
C represent the slope of the magnetic field component in the z and y direction, 
respectively. Inserting (3.9) in (3.5) obtain: 

-f/ S./ H ‘ iyiz = ~^l^l^ f l ‘[ Hi { i ' i+ \' k+x i) +BzJrCv ] ivi! 


= -'‘of (•’■' + 5 ’* + 5)] AzAs < 31 °) 

Note that the magnetic permeability is homogeneous within si since a free-space cell 
was assumed. Its value is thus replaced by the free-space permeability hq. Equating 
(3.8) and (3.10) results in 

- "of [ff,(i,J + i,t + l)]AzAy 

= E, (* . j + 1 . k + 5) - E , (>. j, * + j)] Az 

+ [c,(i, j + i,t) -£„(., i+ifc+l)]As, (3.11) 

Applying central-difference approximation to the time derivative on the left-hand 
side at t = nAt and evaluating the right-hand side of the above equation also at 
t = nAt results in: 



57 


= [£J (i.j + 1 , t + i) - f? (.',;, k + i)] Az 


+ 



(3.12) 


Rearranging and separating new and previous times, provides the equation for up- 
dating the x component of the magnetic "eld in a free-space cell as: 


w * +i ( i,i + 5 ,i + 5) = ifri 0’ i + 5’* + 5) 


+ 


+ 


At 

Ay 

At 

^oA z 


[s" (;, j, * + 1) - e; (», j + 1 . * + i)] 
[£j(i.j + i,i + l)-£;(u + 5.i)] 


(3.13) 


This equation is the same as the one obtained by applying central difference approx- 
imations to the differential form of Maxwell’s equation. 

Applying (3.1) over the Ampere’s contour in Fig. 3.1 results in: 

j> H • dl = ^ J j' fiE x dydz (3.14) 

where Sj is the area enclosed by contour cj. Assuming linear magnetic field com- 
ponents over the four sides of contour C 2 and performing the line integral in (3.14) 


+ [x y (* + \,h k - i) - H v (i + j, k + i)] Ay (3.15) 


results in: 
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Also, assuming linear distribution for the x component of electric field over s 2 and 
performing the surface integral in (3.14) results in: 

J J a tExdydz = ^ e 0 E c (i + ^,j,kJ AyAz (3.16) 

Equating (3.15) and (3.16) yields 

E x (i + ;,*)] AzAy 

= H, (i + i j + 5 ,*) - «. (i + j.J - |.*)] ^ 

+ [i!T» (i + 5 , j, i - i) - /T. (i + 1, >, fc + 1)] Ay (3.17) 

Applying central-difference approximation to the time derivative on the left-hand side 
at (n + 1 12) At and evaluating the right-hand side at (n + 1 / 2 ) A t also, rearranging 
and separating previous and new times results in: 

£,” +1 (;+ j.j',*) = ^(i+j.j.*) 



(3.18) 

The analysis in this section demonstrates that when the integral form of Maxwell’s 
equations is applied over a cubic grid, equations similar to the ones obtained by ap- 
plying central-difference approximations to the differential form of Maxwell’s equa- 
tions are obtained. Since the field components were assumed linear over the cells the 
resulting finite-difference equations are second-order accurate. 

The integral equation formulation has several advantages over the differential 
equation formulation. One advantage is that it provides a better insight as to how the 
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FDTD method works. Also, the integral formulation approach allows the application 
of the FDTD method to a more general class of problems involving distorted grids or 
nonhomogeneous cells. Material discontinuities within the numerical domain can be 
modeled by applying the integral form of Maxwell’s equation at the discontinuities. 

One approach in FDTD applications, is to use modified equations where it is 
only necessary and to apply the standard FDTD algorithm everywhere else within 
the numerical domain. For example, if a curved surface is to be modeled, the grid is 
distorted near the surface of the curved structure only and retains its basic cubic form 
in the rest of the numerical domain. The same approach is applied when modeling sub 
cell features within the numerical domain, for example, very thin dielectric slabs or 
conducting structures coated with thin composite material. To implement correctly 
the required modifications, assumptions concerning the near field components over 
the distorted grid must be made. 

3.3 Contour Path FDTD Method 

Consider now applying the FDTD method to model a smoothly curved conducting 
surface such as the one shown in Fig. 3.2. As illustrated by the figure, the contours 
for updating the magnetic field components must be modified in such a way so as to 
conform with the curved conducting surface. For example, contour Ci is stretched, 
whereas contour C 2 is shrunk to conform with the curved structure. 

Update equations for the magnetic field components within the distorted con- 
tours are obtained by applying Maxwell's integral equations within the distorted 
contours. The electric and magnetic field components are assumed constant within 
the distorted contour segments. This assumption reduces the second-order accurate 
FDTD method to a first-order accurate method. 

The contours for updating the electric fields are not modified. Where possible, 
the electric field components are updated by using normal Ampere’s contours. If that 
is not possible the electric field components are declared unused field points. Par- 




curved conducting 
surface 


Fig. 3.2: Modified contours for implementing the contour path method. 
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ticularly, Ampere’s contours that cross the conducting media interface are avoided. 
For example, the z component of electric field at the left side of contour c 2 cannot 
be updated normally because the magnetic field points are residing on opposite sides 
of the conducting surface. This field component is declared unused field point and 
it is not updated. To update the magnetic field point inside contour c 2 , use is made 
of the lower z component of electric field on the left side of contour c 2 , i.e., a near 
neighborhood approach is used when necessary. 

The application of the contour path FDTD approach to model curved structures 
involves both cubic and distorted grids. The linearity of the field components is 
applied over the cubic part of the grid, while over the distorted part of the grid the 
field components are assumed constant. Therefore, the contour path FDTD approach 
reduces from second-order accurate to first-order accurate over the distorted part of 
the grid. 

The assumptions of the contour path FDTD method can be summarized as fol- 
lows [44], [45]: 

1) The normally rectangular Faraday’s contours surrounding the magnetic field 
components neax a curved conducting surface are distorted so as to conform to the 
surface. 

2) The magnetic field components are assumed to represent the average value of 
the magnetic field within the surface bounded by the distorted contour. 

3) The tangential component of the electric field, E tan , of the distorted contour 
at the conducting curved surface is zero. 

4) Along the straight portions of the distorted contour the electric field compo- 
nents tire assumed constant; where possible these fields are calculated using rectan- 
gular Ampere’s contours from adj acent magnetic field components. 

5) Calculations of Ampere’s contours which cross the boundary of the curved 
conducting surface are avoided; the electric fields along the corresponding Faraday’s 
contour segments, if needed, are obtained using a “near neighbor” approximation. 



62 


The concepts behind the contour path FDTD method will become more clear 
in the next chapter, where the method is applied in the modeling of pyramidal 
horn antennas with or without composite impedance walls. The grid is distorted to 
conform with the flared surface of the horn. The equations updating the magnetic 
fields within the distorted grid are derived based on the assumptions of the contour 
path FDTD method. 

3.4 Contour Path FDTD Modeling of Horn Antennas 

The pyramidal horns were also modeled using the contour path FDTD method. In 
this case the grid was distorted near the antenna surface to conform to its surface. 
The FDTD equations that update the fields near the surface were modified according 
to the distorted grid. These modifications were based on the assumptions stated in 
the previous section on the contour path FDTD method. Away from the antenna 
surface the cubic form of the grid was retained. 

Fig. 3.3 illustrates how the grid can be distorted to model the pyramidal horn 
surface in the El-plane. The FDTD grid is represented using dashed lines. The dots 
in this figure represent the locations of the magnetic field points whose contours are 
affected by the presence of the inclined surface. In the El-plane cross section, in 
which the antenna flares in the y direction, the dots represent x directed magnetic 
field locations. 

Consider the upper surface of the antenna in the El-plane cross section shown in 
Fig. 3.3. For flare angles less than 45°, three distinctive contours were identified. 
These three contours are shown in Fig. 3.4. In Fig. 3.4(a) the antenna surface cuts 
the grid below the location of the magnetic field point in the cell, in Fig. 3.4(b) 
above the location of the magnetic field point, and in Fig. 3.4(c) the antenna surface 
extends from a lower cell point to an upper cell point. 

The equations for updating H x in each case must be modified based on the 
assumptions of the contour path FDTD method. Assume that the electric field 
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Fig. 3.4: Distorted Faraday’s contours used to update H x near the upper surface of the 
horn in the E-plane cross section. 
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components are constant on the sides of the distorted contours and that the magnetic 
field components are constant within the area enclosed by the distorted contour. 
Applying the integral form of Maxwell’s equations over each contour shown in Fig. 
3.4 results in the following update equations for the x component of the magnetic 
field, respectively: 

Hr + * (i,j + k + 2 ) = Hx + 2 ) 

- E; (i,j + i,fc) ■ (I, + Ay) + £" (i,j,k + i) • Az 


(3.19) 


HT 


H iJ 4*4) = ^ + 5. * + 5 ) 


+ 


Ai 

Mo 


E' y 


+ +!)•». 



Az 


(3.20) 






! ?( i,i 4*) 


l2 + E?(i,j,k + ^yAz 


(3.21) 


where A u M, and ^3 are the areas enclosed by the distorted Faraday’s contours (a), 
(b), and (c) of Fig. 3.4, respectively. Note that the contribution from the tangential 
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electric field component on the antenna conducting surface is zero, and therefore, 
not included in the above expressions. 

Similar equations can be written for the distorted contours outside the upper sur- 
face of the horn antenna as well as for the lower surface. Distances l\ and / 2 in each 
case (a), (b), and (c) are shown in Fig. 3.4, and they are estimated from the inter- 
section points of the flared surface with the grid. Once l\ and / 2 are estimated, then 
the areas A \ , A 2 , and A 3 are determined for each case (a), (b) and (c), respectively. 

Assuming a cubic grid, i.e. , Ax = Ay = Az = A, define the following normaliza- 


tions 

n 
1 2 
a; 

A' 


h 

A 

(3.22) 

h 

A 

(3.23) 

A x 

A 2 

(3.24) 

A? 

A 2 

(3.25) 


A' - — 

*3 ~ A2 


(3.26) 


Inserting the above normalizations into (3. 19)-(3.21) simplifies the modified equations 
for updating the x component of magnetic field as follows: 


+ 


At 1 


HoA 


A [£»(<,,- + 1, *+!).(/; + 1.0) 


- e; (*,; + i,t) • (li + 1-0) + e; + 1) 


(3.27) 
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tf " +2 + 2’^ + 2) = 2 (*ii + 2’^ + 2) 


+ 


At 1 


MoA 


A[£ ? (.-,i + i.* + i).r l 


a" + ‘ (.,_/ + j,t + j 


- (<J + J, *) ' G + £? (*,7, * + Jyj 

^•^K( w + 5'* + 0-« + 1J,) 


(3.28) 


+ 


- £j(i,i + i i).r j + £T(M,t + |)] (3.29) 

where the normalized areas for Aj, A^, A 3 enclosed by the distorted contours (a), (b), 
and (c) of Fig. 3.4, respectively, are given by: 


a; = 


(3.30) 

A' = 

1 

2 

(3.31) 

A's = 

1.0 + l[ + 1' 2 

2 

(3.32) 


Finally, the normalized distances l[ and l 2 can easily be determined from the intersec- 
tion points of the antenna surface and the FDTD grid. The above formulation makes 
programming, of the modified equations for updating the magnetic field component, 
easy. 

In certain cases the y directed electric field at the sides of the distorted contours 
c an not be updated because Ampere’s contours that cross the surface of the antenna 
must not be used. The algorithm identifies such cases and uses the most close y 
directed electric field in updating the x directed magnetic field. If the H x field point 
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is located within the antenna surface, the lower E y field is taken. If the H x field 
point is outside and on the upper surface of the antenna, then the upper E y field is 
used. Similar reasoning is applied to the lower surface of the horn antenna. 

For example, consider the distorted grid shown in Fig. 3.5. The electric field 
denoted as E y4 cannot be updated using the magnetic field points H x2 and H x2 , 
because they are at points on opposite sides of the antenna surface and Ampere's 
contours intersecting the boundary media should not be used. Therefore, the value of 
E y4 is not available. To update H x 3 , instead of E y4 use the nearest E y field collinear 
with it; i.e., E y2 . The resulting equation to update H x 3 is 

Hx 3 2 (hi + 2 ' k + 2 ) = 3 (hi + + 2 ) 


+ 


At 1 


Mo 


A[£-( i , i + I,t + l). /l 


- E v 2 (hi + A:) • l 2 + E n zi (i,j, k + ^j-Ai 


(3.33) 

Because the value of E y4 is not available, the z component of the magnetic field 
half-cell above the position of E v 4 cannot be updated. Both the E y4 and H z field 
components are declared unused field points and are not updated by the FDTD 
algorithm. However, the H z magnetic field component is needed to update the field 
points E x 1 and E x2 . These field components are also assumed to have zero value 
(tangential field components close to a conducting surface), and they are not updated 
by the FDTD algorithm. 

Other instances where some of the distorted Faraday’s contours create field points 
that should neither be advanced nor be used by the algorithm to advance other 
neighboring field points are now examined. Examples of unused field points can 
easily be identified in case (c) of Fig. 3.4. In this case the field points corresponding 
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Fig. 3.5: Example of modified updating at locations where an Ampere’s contour cannot 
be used because it is crossing media boundaries. 
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to the y-directed magnetic field and the z-directed electric field, where the inclined 
surface intersects the upper side of the unit cell, must be declared unused field points. 

The corresponding modifications in the H-plane cross section of the pyramidal 
horn are obtained by following a similar procedure as that applied in the E-plane cross 
section. Fig. 3.6 illustrates how the grid can be distorted to model the pyramidal 
horn surface in the H-plane. The dots in this case represent the locations of the y 
directed magnetic field components whose contours must be modified to conform to 
the antenna surface. The equations updating the y-directed magnetic field points 
in the distorted grid are obtained by applying Maxwell’s integral equations over the 
distorted grid. The procedure is similar to that applied for the E-plane cross section 
and will not be discussed further. 

Slightly different modifications are required for the z component of the magnetic 
field, compared to those applied for the x and y components of the magnetic field. 
Fig. 3.7 illustrates how the grid can be distorted to model the pyramidal horn in the 
xy plane. The dots in this figure represent the locations of the magnetic field points 
whose contours are affected by the presence of the antenna conducting surface. 

Consider the upper surface of the antenna shown in Fig. 3.7. There are two 
possible ways that the antenna surface can cut the grid. First, the antenna surface 
can cut the grid below the magnetic field point and second, above the magnetic field 
point. These two cases ace illustrated in Fig. 3.8. The equations for updating H z in 
each case must be modified by applying the contour path method over the distorted 
contours. Normalizing the distance l\, shown in Fig. 3.8, with respect to the cell 
size and applying the integral form of Maxwell’s second equation (3.2) over the two 
distorted contours, respectively, results in: 

H * +h i i + \ J + \' k ) = + + 
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Fig. 3,7: Distorted Faraday’s contours at the antenna surface in the xy plane. 
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where the normalized areas A[ and A' 2 enclosed by the distorted contours C\ and ci 
of Fig. 3.8, respectively, are given by: 


+ 

0 

II 

(3.36) 

A' - /' 

/± 2 — i x 

(3.37) 


The normalized distances l[ and 1' 2 can easily be determined from the intersection 
points of the antenna surface and the FDTD grid. 

3.5 Input Power, Radiated Power and Radiation Efficiency 

In analyzing the pyramidal horn, it is assumed that it is fed by a rectangular waveg- 
uide operating in the dominant TE\q mode whose tangential electric field at the 
aperture is represented by 

E v (x,y,z) = E 0 sin(-x)sin(ut — 0 x z) (3.38) 

a 

where 0 < x < a and 0 < y < b. The constant /3 X represents the waveguide 
propagation constant. At the reference feed plane j3 x z was set equal to zero. 
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Using this feeding scheme, the input power to the horn was estimated by inte- 
grating the input power density over the waveguide cross section, and it is givenn 

^ (£) (t) - <**> 

where r? 0 = 120x is the free- space impedance. 

The far zone Eg and E$ electric fields are obtained from the FDTD code through 
a near-to-far-field transformation. The El-plane gain pattern was calculated using 

f%(M=90*) 


G e (0, <j> = 90°) = 4x- 




ptnput 


(3.40) 


Similarly, the H-plane gain pattern was calculated using 


EQ»,+= o*) 


(3.41) 


G«(M = 0°) = 4x-^T 

The antenna radiated power was estimated by integrating the far zone fields over 
a sphere and it is represented by [4]: 

t a* 


'--rf U(0, <j>) sin 0d9d<f> 


(3.42) 


where U ( 0 , <j > ) is the intensity given by 


V(!A) = T-ll E'(»A) I* + I l J ] 

ZT)o 


(3.43) 


One important design parameter of pyramidal horns with lossy materials on the 
E-plane walls is the power loss or the antenna efficiency calculated using 

prad 


V = 


pinpvt 


(3.44) 


3.6 Numerical Results 

The 20-dB standard gain pyramidal horn antenna was also analyzed using the con- 
tour path FDTD. Figs. 3.9 and 3.10 compare the computed E- and H-plane gain 
patterns, respectively, of the pyramidal horn with measurements. As illustrated by 



76 


both figures, there is very good agreement between the computed and measured re- 
sults over a dynamic range of 60 dB. The computed gain patterns agree very well 
with the measured ones when a standard cell size of A/12 was used (0.1"). The 
overall FDTD grid required for this problem was 76 x 60 x 142 cells. This simula- 
tion was run for 40 cycles to reach steady state and took about 2,500 seconds on an 
IBM-RISC/6000 mainframe computer. 

The antenna radiation efficiency was estimated by integrating the far- zone fields 
to obtain the radiated power. The amplitude of the y-directed electric field was 
set equal to Eq = 1 V/m, producing a reference power of P\q = 0.116226 x 10 -6 
Watts. By estimating the antenna far fields and then integrating them (1° steps) 
the radiated power was estimated P Ta< i = 0.115077 x 10 -6 Watts. The corresponding 
radiation efficiency was rj = 0.99011 ~ 99.0%. Therefore, most of the input power is 
radiated except for a small portion which is reflected back into the antenna. 

To examine further the accuracy of the contour path FDTD method, it was 
also used in the analysis of the two square aperture pyramided horn antennas. The 
aperture dimensions of the smaller square horn were 5" x 5", whereas the aperture 
dimensions of the larger one was 7" x 7". 

In Figs. 3.11 and 3.12 the computed antenna gain patterns are compared with 
measurements for the smaller square aperture horn. As illustrated by both figures, 
there is very good agreement between the computed and measured antenna gain 
patterns. The frequency of operation was 10 GHz. The aperture of the horn was 
5" x 5" whereas the waveguide aperture was of standard X-band dimensions; i.e., 
0.9" x 0.4". The transition length from the waveguide aperture to the horn aperture, 
L, was 10.5". A grid size of 0.1" was used, i.e., the waveguide was 9x4 cells whereas 
the horn aperture was 50 x 50 cells. The distance from the waveguide transition to 
the horn aperture was 105 cells. 

The computer code was run for 40 cycles to reach steady state, on an IBM- 
RISC/6000 computer. The CPU time was 4,200 seconds for obtaining the El- and 
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Fig. 3.9: E-plane gain of a 20-dB pyramidal horn at 10.0 GHz (A = 4.87", B 
= 10.06", a = 0.9" and b = 0.4"). 
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Fig. 3.10: 
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H-plane patterns over 360° at 1° steps. An overall grid size of 78 x 78 x 144 cells was 
used for this simulation. 



Fig. 3.11: E-plane gain of a square aperture pyramidal horn at 10.0 GHz (A = 5", B = 5", 
L = 10.5", a = 0.9" and b = 0.4"). 

Finally, the contour path FDTD method was used to model a larger square aper- 
ture antenna. The aperture of the horn was 7" x 7" whereas the waveguide aperture 
was of standard X-band dimensions; i.e., 0.9" x 0.4". The transition length from the 
waveguide aperture to the horn aperture, L, was 12.2". The frequency of operation 
was 10 GHz. A grid size of 0.1" was used; i.e., the waveguide was 9x4 cells whereas 
the horn aperture was 70 x 70 cells. The distance from the waveguide transition 
to the horn aperture was 122 cells. Figs. 3.13 and 3.14 compare the computed E- 
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and H-plane patterns of the square aperture horn antenna with measurements. As 
illustrated in the figures, there is very good agreement between the computations 
and the measurements. 
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Fig. 3.13: E-plane gain of a square aperture pyramidal horn at 10.0 GHz (A = 7", B = 7", 
L = 12.2", a = 0.9" and b = 0.4"). 








CHAPTER 4 


CONTOUR PATH FDTD ANALYSIS OF HORN ANTENNAS WITH LOSSY 

WALLS 


4.1 Introduction 

An advantage in the use of the FDTD method is its ability to model thin composite 
materials such as coated conductors with dielectric and/or magnetic material slabs. 
Such materials find applications in electrically and/or magnetically coated aperture 
antennas [40]. 

Analysis of three-dimensional material plates, with and without conductor back- 
ing, was presented using an open surface integral formulation by Newman and Schrote 
[46]. Min et al. [47] presented an efficient formulation, using the integral equation ap- 
proach, to solve scattering problems from two-dimensional conducting bodies coated 
with magnetic materials. A hybrid finite element method has been applied to con- 
ducting cylinders coated with inhomogeneous dielectric and/or magnetic materials 
[48], and a combined finite element-boundary integral formulation was developed for 
the solution of two-dimensional scattering problems coated with lossy material [49]. 

The thin coating of composite material presents a computational problem when 
the FDTD method is applied, because the method demands that the grid size used 
is at least equal to the thickness of the coating. When the grid size of the FDTD 
method is chosen very small, the corresponding time step increment must also be 
comparatively small to satisfy the stability condition [50]. This requirement results 
in tin inefficient way to solve the problem. 

A modified FDTD formulation based on the integral representation of Maxwell’s 
equations wa s applied to thin dielectric structures [51]. The structures involved 
thin dielectric slabs (with and without conductor backing) and conductor- backed 
dielectric slabs that contain thin dielectric cracks in the coating. The dielectric 
coatings considered were lossless, i.e., the conductivity of the dielectric material was 
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zero. Using a slightly different approach, Maloney et a 1. [52] applied the FDTD 
method to thin dielectric slabs with finite conductivity. 

In this research project, the formulations presented in [51] and [52] were extended 
to model pyramidal horns coated either partially or totally by a thin, lossy, magnetic 
material coating. The magnetic coating, having both electric and magnetic loss, 
was applied in the E-plane upper and lower walls of the antenna for pattern control 
purposes. Depending on the extent to which the E-plane wall is coated, nearly 
symmetric E-plane and H-plane antenna patterns can be obtained. The approach 
followed to include the effect of such lossy materials is described next. 

4.2 Contour Path FDTD Modeling of Horn Antennas With Lossy Walls 

The geometry of a partially coated pyramidal horn antenna is shown in Fig. 4.1, along 
with the E- and H-plane cross sections of the antenna. As illustrated by the figure, 
the upper and lower walls of the horn axe coated on the inside with a thin composite 
material having both electric and magnetic losses. Since the composite material 
is placed only on the upper and lower E-plane walls, modifications in the FDTD 
update equations must be made only in the E-plane cross section of the antenna. 
Also, since the composites are placed only on the inside part of the pyramidal horn, 
the equations for the inside part only, need further modifications (from the metallic 
case). The modifications made for the outside part of the metallic horn apply for 
this case also. 

Fig. 4.2 shows how the presence of a thin composite material on the inside of 
the upper and lower El-plane waills modifies the FDTD grid. The composite material 
section is shown as the shaded region. This region covers, in some cells, the magnetic 
field points, shown as dots, whose contours are affected by the presence of both 
the antenna surface and the lossy material. Like in the metallic horn case, three 
distinctive contours were identified. These contours, for the upper section of the 
wall, including the lossy coating, are shown in Fig. 4.3. 
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Fig. 4.1: Geometry of partially coated pyramidal horn antenna. 
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Fig. 4.3: Distorted Faraday’s contours used to update H x in the presence of a thin section 
of composite material. 

The distorted contours shown in Fig. 4.3 are partly-filled with lossy material 
and partly-filled with free space. The derivation of the modified equations for this 
case will involve applying Maxwell's integral equations within the distorted contours. 
Since the contours are partly-filled with material and partly-filled with free space, 
care must be taken to ensure that the physics of the problem is not violated. 

Consider first deriving the equation for updating the magnetic field component 
in case (a) of Fig. 4.3. Applying (3.2) over the distorted contour results in: 

i Edl = -I//, (41) 

where C\ is the contour shown in Fig. 4.3(a) and Si is the total area enclosed by the 
contour C\ . Splitting the integrations over the material part and free-space part of 
the cell and taking the partial derivative with respect to time within the integral, 
the following equation is obtained: 

<f E-d\ = - [ [ da - [ f pH da 

Jci J 1-imj Ot J 
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where Si — s mi denotes the free-space part of the cell and s m , the material part. The 
next step is to perform the integrations in (4.2) based on the assumptions of the 
contour path method. 

First, consider normalizing distances /i, I 2 and d with respect to the cell size A 
and the total distorted area si and material area s mi with respect to the square cell 
area A 2 . This normalization results into corresponding primed variables l[, 1' 2 , d! , s\ 
and s' m . Using these normalizations, the total normalized area of the cell is given 
by: 

/' 4 - 1' 

s\ = 1 4- (4.3) 

whereas, the composite material area is given by: 

*:», = * (4.4) 

Assuming constant magnetic field distribution for H x within the distorted cell 
and integrating the right-hand side of (4.2) results in: 

E - tfl = - W fl + ^ - /) A 1 ^ - pcP,<t - p<r^H, 

= + (4.5) 

where \t T and p are the relative dielectric permeability and magnetic resistivity, re- 
spectively, of the composite material. Introducing the parameters 

TT, = 1 + - <r + Pr* (4.6) 

f = pi (4-7) 

simplifies (4.5) into the following: 

E • <A = -pojSA 1 ^ - ?A >H, (4.8) 


O' 2. 
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To perform the line integral of the electric field on the left-hand side of the above 
equation, assumptions must be made about the distribution of the y component 
of electric field in the free-space part and material part of the distorted contour. 
Since E y is not a tangential field, at the air/material interface its value undergoes a 
discontinuity. This discontinuity must be included in the model. 

The assumed distribution of the y component of electric field on the right-side 
of contour Ci is shown in Fig. 4.4. The discontinuity of the field at the air/material 
interface is shown in this figure. E y i n represents the value of the field within the 
material and E yont the corresponding value in the free-space part of the cell. 



Fig. 4.4: Assumed distribution for E y on the right side of contour c\. 

To perform the line integral of the y component of electric field over the right 
side of the cell, E yxn and E yout must be integrated over the material section and 
free-space section of the right side of contour ci, respectively. A similar treatment 
must be made for the left side of the contour c\. The line integral of the electric field 
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around contour c\ can be expressed as: 


/ 
Jc 1 


E • d\ 


I E am , (1.0 + I'l - i) A + £* n <fA]„ fl 


- [«, - (i.o + r 1 -fl 4 + e n ,MU-M| lm (4.9) 


This formulation requires two values for the y-directed electric field component, 
one in the free-space part of the distorted cell and one in the material part. The 
FDTD algorithm however, updates only one value for the y-directed electric field on 
each side of the distorted contour. This updated field component is either assumed 
in free space and hence the free space parameters are used in updating it or in the 
material and hence the material parameters are used in updating it. To alleviate 
the problem of having to estimate two electric field values for the y component of 
electric field in the distorted part of the contour, the value of the electric field inside 
the material is related by the corresponding value in free space using a boundary 
condition relation at the air/material interface. 

The additional expression required to relate E v ,„ and E you t is obtained by applying 
the boundary condition for E yin and E yout at the air/material interface. Note that in 
applying the boundary condition, the inclination of the interface must be accounted 
for. The additional expression from applying the boundary condition for E y at the 
air /material interface can shown to be: 


e l sin 2 0 e + cos 2 0 e 


F — 

-‘-'yin — 


E 


tr 


yoxit 


(4.10) 


where t r is the relative dielectric permittivity of the material and 6 e is the E-plane 
flare angle of the antenna. 

Assuming that both the y-directed electric field components in contour c\ are 
updated in free-space medium, the line integral of the electric field is evaluated by 
replacing the value of Eyin in terms of the E you t in (4.9). This substitution results in 
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the following expression for the line integral: 


f ci - {Ev l 2‘ff A ] Uft [ E v l 'ufA} rtght - [E,A] bottom (4.11) 

where E yout was replaced by E y and l u!f and 1'iejj are given by the following expres- 


d' del sin 2 0 r + cos 2 # e 

I'm, = (1.0 + /;-/) + -V-: — j : 


(4.12) 


/ , , ^ V sin 2 + cos2 # 

he}} = (1-0 + / 2 — (f ) H 


(4.13) 


Replacing the derived expression for the line integral (4.11) into (4.8) results in the 
following simplified expression: 

a] wi - [£y; t „A] r . jti - [£,A] k , lOT , = -^a^ - - p cSh x (4.h) 

The above equation can now be discretized by following the standard procedure 
of the FDTD method. The final expression for updating the x component of the 
magnetic field in the distorted contour of Fig. 4.3(a) is: 

i pAt 

+ !,* + ») = - -W + + 

r I* n . — — 


where 


A t_ 1 
MOMr A 

1 + 


[f? (■•> + J.* + 0 I'm, - K ('.2 + k ) 

«?(«•*+ 5 )] 


(4.15) 


1 + 4+4_/ +m - 


(4.16) 
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= pd' 


(4.17) 


^le / / ~ (1-0 + l[~(?) + 


d! sin 2 9 e + cos 2 0 e 


(4.18) 


l 2eff ~ l 1 ' 0 + l 2 ~ d ') + 


d'\J e* sin 2 9 e + cos 2 9 t 


(4.19) 


Following a similar procedure, equations for updating the x component of the 
magnetic field in contours c 2 and c 3 of Fig. 4.3(b) and (c), respectively, can be 
derived. For contours c 2 and c 3 , (4.15) applies also, but with different effective 
electrical parameters and effective lengths. 

The expressions of effective electrical parameters and lengths for updating the x 
component of the magnetic field in the distorted contour of Fig. 4.3(b) are given by: 


/' -l- V 

Tr = ^ ~ d! + tird! 


(4.20) 


= pd! 


(4.21) 


'w/ = ('1--0 + 


d'\fd* sin 2 0 t + cos 2 6 e 

€ r 


(4.22) 


tu, = («,-*) + 


d'yjtl sin 2 9 e + COS 2 9 e 

(■T 


(4.23) 


The expressions for updating the x component of the magnetic field in the dis- 
torted contour of Fig. 4.3(c) are: 

1.0 + 1 \ + f 2 j, - /, 

T = 7T - dt + p T dC (4.24) 


= pi 


(4.25) 


I'un = (1.0 + r,-<0 + 


d'yjt* sin 2 0 C + cos 2 9 t 


(4.26) 
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*W/ = ( l 2 -<*') + 


d! J e 2 sin 2 9 e + cos 2 9 e 


(4.27) 


Slightly different modifications are required for the z component of the magnetic 
field in the E-plane cross section, compared to those applied for the x component of 
the magnetic field. Fig. 4.5 illustrates how the grid can be distorted to model the 
pyramidal horn in the xy plane. The dots in this figure represent the locations of 
the magnetic field points whose contours are affected by the presence of the antenna 
conducting surface and composite material. In this case, the distorted contours are 
normal to the air/material interface. Thus, the boundary condition for E yyn and 
Eyout at the air/material interface can be obtained by setting 9 e = 0° in (4.10). The 
equations for updating H z are then obtained by following a similar procedure as the 
one applied to the x component of magnetic field. 

Consider the upper surface of the antenna shown in Fig. 4.5. There are two 
possible ways that the antenna surface can cut the grid. First, the antenna surface 
can cut the grid below the magnetic field point and second, above the magnetic field 
point. These two cases are illustrated in Fig. 4.6. The equations for updating H z in 
each case must be modified by applying the contour path method over the distorted 
contours. Normalizing the distance l\, shown in Fig. 4.6, with respect to the cell size 
and applying (3.2) over the distorted contour Ci results in: 

. 1 _ , 
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MO Mr A 

1 + 

1 T 


[f? («.j + *) - e; (i + ij + i, t) i\„ 


- £;(i + i.>.t)] 


(4.28) 




Fig. 4.5: Distorted Faraday’s contours at the antenna surface in the xy plane with the 
presence of composite material. 
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Fig. 4.6: Distorted Faraday’s contours used to update H z near the upper surface of the 
horn in the xy plane with the presence of composite material. 
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where 


/i r — 1 + — d! + fi r d? 

(4.29) 

P = P<? 

(4.30) 

= (i +/;-<o + 7 

(4.31) 

For contour c 2 of Fig. 4.6(b), (4.28) also applies. The effective electrical parame- 
ters and effective length in this case are slightly different. The parameters for contour 

C 2 are 


/j, r — /j — df -f - /ir(f 

(4.32) 

p = p<? 

(4.33) 

l'u„ = «-<0 + 7 

(4.34) 


The modifications for H x and H x were derived based on the assumption that 
the material thickness was much smaller than the distorted FDTD cells near the 
antenna surface. By increasing the material thickness but still keeping it less than 
the FDTD cell some slight modifications must be made. For example, consider the 
three distorted contour shown in Fig. 4.7. The evaluation of the electric field line 
integral in this case depends on whether the y-directed electric field is updated using 
free- space parameters or composite material parameters. If the FDTD updated value 
for E y used material parameters the value of E vou t in (4.9) is substituted by E yin 
using the air/material boundary condition (4.10). 

In Fig. 4.7(a) the y-directed electric field components on the left and right sides 
of contour ci are updated using free-space parameters. Hence the expressions derived 
in the previous case with thin material also apply in this case. For Fig. 4.7(b) both 
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(a) (b) (c) 


Fig. 4.7: Distorted Faraday’s contours used to update H x in the presence of a thick section 
of composite material. 


y-directed electric field components are updated using material parameters. Hence, 
the following expressions for the effective electrical parameters and lengths are used: 


Mr = 

/' 4- /' 

^P-^ + Mr^ 

(4.35) 

P = 

p<r 

(4.36) 

/' _ 
‘w/ - 

7 (r - +<r 

y ( cj sin 2 6 e + cos 2 9 e 

(4.37) 

hejf ~ 

& *+<t 

\J e 2 sin 2 0 e + cos 2 $ e 

(4.38) 

For Fig. 4.7(c) the y-directed electric field component on the left side of contour 
C 3 is updated using material parameters, and the corresponding commponent on 
the right side using free-space parameters. Hence, the following expressions for the 
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effective electrical parameters and lengths are used: 


/' -1- /' 

Tr = 1 2 2 d' + H T d! 

(4.39) 

p = pdf 

(4.40) 

dfJ e 2 sin 2 9 e + cos 2 9 e 

^ = (i.o + /;-<o + — — ~r 

(4.41) 

r 2ef{ - (/'-(f)- 7= — Cr +<f 

\j e 2 sin 2 9 e + cos 2 9 e 

(4.42) 


The above modifications constitute the basis of modeling pyramidal horns with 
composite E-plane walls using the contour path FDTD method. The assumptions 
made to derive these modified equations were based on the contour path FDTD 
method. For making possible such am implementation, certain assumptions concern- 
ing the field distribution of electric and magnetic fields in the distorted grid were 
made. Because of these assumptions, the above algorithm would produce accurate 
results under certain conditions. The limitations of the contour path FDTD method, 
when applied in the analysis of pyramidal horn antennas with and without composite 
E- plane walls, is the subject of the following section. 

4.3 Limitations of the Contour Path FDTD Method for the Analysis of Radiation 
by Pyramidal Horns 

The flared surface of pyramidal horns was modeled by distorting the FDTD grid near 
the antenna surface. The integral form of Maxwell’s equations was applied over the 
distorted grid to obtain modified update equations for the magnetic field components 
within the distorted grid. The field distribution over the distorted grid was assumed 
constant. This assumption reduces the FDTD method from a second-order accurate 
method in space to a first-order accurate method. The method, however, remains 
second-order accurate in time. 
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The decrease of one order of accuracy in space is expected to have some effect 
on the accuracy of the computations. The reduced accuracy is expected to affect 
more the modeling of pyramidal horns with composite E-plane walls. Because of the 
presence of the composite material (with high dielectric permittivity and magnetic 
permeability constants and high magnetic loss), the field distribution is expected to 
vary significantly within the material. However, when the thickness of the composite 
material is electrically small (less than A/20), the constant field assumption over 
the cell will be more accurate because of the presence of the conducting surface. 
The slope of normal electric field components, for example, is zero near a perfectly 
conducting surface. 

The implemented contour path FDTD method is valid for the analysis of com- 
posite materials with thickness smaller than the FDTD cell. Should the composite 
material thickness become larger than a unit cell, a different implementation would 
be required. 

Another assumption made in implementing the contour path approach was that 
the antenna flare angles in the E-plane and H-plane cross sections, were smaller 
than 45°. With flare angles larger than 45° a more complex algorithm has to be 
implemented. The three distinctive contours in Fig. 3.4 were identified based on the 
assumption that the flare angles were less than 45°. More possible contours could 
be identified, with the angle being larger than 45°. Hence, the algorithm complexity 
would increase. 

4.4 Numerical Results 

In this section the developed FDTD computer program is applied for the analysis 
of pyramidal horn antennas with composite E-plane inner walls. The numerical 
results are validated by comparing with available experimental data. The effect of 
composite material on the gain pattern of the antenna is examined and illustrated 
through a number of computed patterns. Radiation efficiencies and design curves 
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for the variation of the broadside gain loss as a function of material thickness and 
length are also presented. 

The standard 20-dB gain horn antenna is analyzed first. For pattern control 
purposes sections of ECCOSORB GDS composite material with measured electrical 
parameters e r = 14.9 — j0.25 and /i r = 1.55 — jl.45 at 10 GHz were placed on the 
inner E-plane walls of the horn as illustrated in Fig. 4.1. The nominal thickness 
of the composite material section was 30 mil (0.0762cm). The material thickness 
however, was measured to be t = 33 mil (0.08382cm). Two different material lengths 
l m were used for this case. The first case corresponds to l m = 2" (5.08cm) and the 
second case to l m = 4” (10.16cm). 

For the FDTD simulations a grid size of 0.1" was used. The material thickness 
was about a third of the FDTD grid. The composite material sections influence 
mainly the El-plane pattern. Thus, only results from the El-plane gain pattern calcu- 
lation of the coated 20-dB standard gain horn are presented. Fig. 4.8 compares the 
FDTD computed El-plane gain of the antenna with measurements, when sections of 
2" composite material are used. The agreement between computed and measured 
results is good. As illustrated in the figure the first sidelobe of the pattern is ellimi- 
nated. Because of the presence of magnetic material in the inner walls, the broadside 
antenna gain was reduced by 2.69 dB. The antenna radiation efficiency for this horn 
was found to be 77 = 0.76315 — 76.3%. 

The composite material is more effective in terms of reducing the antenna side 
lobe levels, when a larger section of composite material is used. Fig. 4.9 compares 
the FDTD computed results with measurements , when a 4" section of composite 
material was used. The agreement between the computed and measured results is 
good. As illustrated in the figure, the second side lobe is now almost eliminated. The 
side lobe levels in the back side of the antenna are reduced to levels below -10 dB. 
The reduction in the broadside antenna gain in this case was close to 5.0 dB. The 
radiation efficiency was reduced to 77 = 0.63870 ^ 63.9%. In this case a significant 
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amount of power is dissipated in the lossy magnetic material. 
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Fig. 4.8: El-plane gain of a 20-dB standard gain pyramidal horn at 10.0 GHz, partially 
coated with GDS magnetic material (e r = 14.9 — j'0.25 and \i T = 1.55 — j 1.45, t = 33 
mil and l m = 2"). 


To examine the effect of material thickness and length on the antenna radiation 
pattern, the broadside gain loss of the 20-dB standard gain horn was calculated 
for different GDS material thicknesses and lengths. The variation of the broadside 
antenna gain loss versus the GDS material thickness, for different material lengths, is 
illustrated in Fig. 4.10. As illustrated in the figure, the broadside antenna gain loss 
is largest for thickness in the range of 25 - 30 mil and decreases for larger material 
thickness. For the larger thickness range the gain loss is not strongly influenced by 
the material length. 




102 



Fig. 4.9: El-plane gain of a 20-dB standard gain pyramidal horn at 10.0 GHz, partially 
coated with GDS magnetic material (e r = 14.9 — j0.25 and p r = 1.55 - j‘1.45, t = 33 
mil and l m = 4"). 
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The effectiveness of the thicker material was demonstrated by two examples. Two 
sections of 33 mil GDS material were combined to form a 66 mil thickness. The 66 
mil thick material was then applied to the 20-dB standard gain horn antenna. Two 
different lengths of material were used in this case. In the first case the material 
length was 2" and in the second case 4". 

The computed E-plane gain pattern of the 20-dB standard gain horn antenna, 
with a 2" long and 66 mil thick GDS material on the inner E-plane walls, is com- 
pared with measurements in Fig. 4.11. As illustrated in the figure, there are some 
differences between the computed and measured antenna gain patterns. These dif- 
ferences are attributed to the larger thickness of the composite material. With a 
thicker section of composite material on the antenna inner E-plane walls, the as- 
sumption of constant field distribution in the vicinity of the distorted contours is no 
longer accurate. However, both the FDTD computed results and the measurements 
demonstrate in this case that the 66 mil thick material case is more effective in shap- 
ing the antenna pattern, with out considerable broadside gain loss and reduction in 
the antenna efficiency. The broadside gain loss for this case was about 2.6 dB, and 
the antenna efficiency 77 = 0.70 70.0%. 

The results of the second case with 4" length and 66 mil GDS material are exhib- 
ited in Fig. 4.12. In this case also there are differences between the computed and 
measured antenna gain patterns. The reduction in the broadside gain for this case 
was about 2.8 dB, and the antenna efficiency was tj = 0.78559 ~ 78.56%. 

4.5 Conclusions 

The contour path FDTD method was used to analyze pyramidal horns with or with- 
out composite E-plane walls. For the metallic horns accurate gain patterns were 
computed and compared to available measurements. When the inner El-plane walls 
were coated with lossy materials, for gain pattern control, the contour path FDTD 
method yields acceptable gain patterns compared to measured ones. For this imple- 
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Fig. 4.10: Broadside antenna gain I 068 of a 20-dB standard gain pyramidal horn at 10.0 
GHz, partially coated with GDS magnetic material (c r = 14.9 — j0.25 and fi T = 1.55 — 
jl.45). 
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Fig. 4.11: E-plane gain of a 20-dB standard gain pyramidal horn at 10.0 GHz, partially 
coated with GDS magnetic material (f r = 14.9 - j'0.25 and /i r = 1.55 — jl.45, t = 66 
mil and l m = 2"). 
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Fig. 4.12: E-plane gain of a 20-dB standard gain pyramidal horn at 10.0 GHz, partially 
coated with GDS magnetic material (f r = 14.9 - j0.25 and fi r = 1.55 - j'1.45, t = 66 
mil and l m = 4"). 
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mentation some limitations were imposed. 

One approach to design horn antennas with low side lobes and rotationally sym- 
metric patterns is by using lossy materials. This approach was investigated during 
this project using the FDTD method. Depending on the electrical properties, length 
and thickness of the material used, nearly symmetric antenna patterns can be de- 
signed. Due to the presence of the lossy material, the broadside gain of the horn 
is reduced by few dB’s. However this loss of gain can be minimized by properly 
selecting the material dimensions, especially the thickness. Design of antenna gain 
patterns based on the developed computer program is now possible. Such computer- 
aided designs produce a desired pattern with minimum reduction in the antenna 
broadside gain. 
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APPENDIX A 

EVALUATION OF THE IMPEDANCE MATRIX ELEMENTS IN MOMENT 

METHOD 
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The electric fields due to surface electric and magnetic current distribution needed 
in (1.48) and (1.49) can be written as 


E(J) = -jfivj/ + (A.l) 

E(M) = - / M(«') x V$ds' (A. 2) 

J s* 

where primed coordinate represents the source coordinate and $ = e ^ R / 47t R is 
the free-space Green’s function; J(s') and M(s') are surface electric and magnetic 
current density on S', respectively. Substituting (A.l) and (A. 2) into (1.48) and 
(1.49), respectively, using the current continuity condition, the impedance matrix 
elements can be expressed in the form of 

= 113 J s L l{p ‘ ' PjJ) - ^ (V • P ' ){V ■ p/)]**;*** (*•») 

<?« = /»// ^PiMPFM )(1 +jf>R)i<ts',d st (A-4) 

Js k JS' t nr 

where n is the unit directional vector of (P^ x Pj*). Evaluations of (A.3) and (A. 4) 
can be efficiently carried out by using Taylor’s expansions of the kernel functions $ 
and e~ 10R at the center of the two patches and analytical integrations can be found 
for individual terms of the Taylor’s expansions. Therefore, only one of the surface 
integrals needs to be evaluated numerically. 







